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ABSTRACT 


A seven level primitive equation equatorial ocean model 
with temperature and salinity is integrated over a 30 year 
period to a quasi-steady state. Placement of four of the 
seven levels of the model in the first 100 meters gives 
sufficient vertical resolution to resolve most major features 
of equatorial oceans. The model equivalents of the South 
Equatorial Current, the Cromwell Current, vertical upwelling 
in the equatorial area, Ekman divergence of mass toward each 
s and aerriermee line were all simulated. 

A major conclusion was salinity probably has little 
effect on the general circulation patterns in the equatorial 
region. Additionally, the use of a vertical eddy diffusion 
coefficient for momentum that is dependent upon vertical 
stability, vertical current shear, and depth will probably 
be required to accurately model the vertical structure of 


ehe currents. 
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|. INTRODUCTION 


Numerical baroclinic ocean models are advancing to the 
state where oceanic forecasting may become feasible. The 
motivation for the development of sophisticated ocean models 
sbems largely trom the meed Co aceurately speciiy since vc. 
Boundary conditions for baroclinic atmospheric models used 
in long-range weather prediction. The importance of the 
ocean in determining global climate is underscored in a 
recent pioneering work by Manabe (1969) and Bryan (1969) in 
which numerical ocean and atmospheric models were coupled 
together allowing interaction between the two. By comparing 
the results of the ocean-air model with those from an atmos- 
pheric model with a static ocean, they found the ocean plays 
an important role in the poleward transport of heat and 
-helps determine the thermal structure and vertical wind shear 
of the troposphere. Additionally, they found that the ocean 
affects the hydrology of the tropics and influences the 
eereneth of the Hadley circulation. The important role of 
the ocean in influencing the transient behavior of the 
atmosphere has been recently shown by Lambertson (1972) in 
a coupled ocean-atmospheric model. Because of the demon- 
arated importance ol the ocesu ти Ссил ии ии ии e 
and the variability of the atmosphere and sia eramos pcs 
models are more highly developed than their ocean counter- 
parts, it is clear that future improvement in long-range 
weather prediction is closely linked with the development of 


better ocean models. 





Barly oceanic prediction moael5 were one mona iona 
the vertical with the net heat flux and, in some cases, the 
wind velocity or stress being specified at the surface. 

The thermal structure was assumed to be known initially 

and responded to changes in the rates of surface heating. 
The depth of the surface TNES was determined by a mechan- 
ical energy balance, and integration of the thermodynamic 
equation gave the thermal profile at a later time. Although 
adequate for very short range purposes, these one-dimensional 
Models are not satisfactory for long range predictions 
because no mechanisms for horizontal advection or diffusion 
of heat or momentum are present. Any location where the 
climatological yearly heat flux is not zero will display 
either a continual warming or cooling trend. Tropical 

areas would continue to warm ana middle latitude areas 

oc continue Co cool. Por long term predictione 201657 
Mare. ocean models based on heat balances alone will not 

Be satisfactory. 

Although Bryan's dynamical ocean model is adequate for 
mae purpose of global climate investigation, it 15 believed 
that the horizontal grid spacing of 5° longitude and 5° 
latitude and a vertical spacing where the first level is 
at 67 meters is too coarse for adequate oceanic forecasting 
purposes. One strong point of the present model is felt 
to be its superior horizontal resomutten ver 2.5° Lav cud. 
and 3.0? longitude and good vertical resolution where four 


of its seven levels are in the first 100 meters. 
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Because of the importance of equatorial air-sea inter- 
action (Bjerknes, 1966) and the interesting oceanographic 
phenomena that occur in the equatorial region, it was 
decided to model the equatorial domain. The large vertical 
current shears and sharp thermocline which exist in the 
tropics should be simulated by an ocean model that is to 


be used with an atmospheric forecasting model. 


TE 





IT, DESIGN OF THE MODER 


A. MATHEMATICAL STATEMENT OF PROBLEM 

The idealized domain is 60° wide, extends from 21.25 S 
вә 21.25 N, and has a uniform depth of 358000 meters.  Plaece= 
ment of the lateral walls at approximately 20° latitude 
in both hemispheres satisfies two important requirements: 
first, the curl of the observed wind stress is almost zero 
ae chese locations: and Second, Che arca averageqe: luxvor 
salt as determined from climatological evaporation and 
Poecipitatiom data is approximately zero when averaged Trom 
21.25 S to 21.25 N. Seven vertical levels placed at 5, 20, 
"or 100, 250; /00, and 2000 meters are used. The upper 
levels are more closely spaced to provide better resolution 
=n Che near surface domain where the largest variation in 
dependent variables is found. Horizontally the grid points 
are spaced 3.0 degrees longitude and 2.5 EEE latitude. 
ime horizontal and vertical domains are pictured in figures 
51 and 2.2 respectively. 

This model is similar to the primitive equation models 
of Haney (1971), Bryan (1969), and Bryan and Cox (1968) 
except forthe introduction of salinity ср ии ас 
pendent variable. The solution of the equations of motion 
involves the use of a streamfuncilon 1705756507 0. 0) 
motion. This somewhat artificial constraint a Ove Ome ии 


of a large time step and saved considerable computation time. 
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Figure 2.1 Horizontal domain of the model Horizontal 
velocity components m and v were carried between grid 
points. The variables Y, 5, », ¥,w Were carried on tae 
grid points. Grid spacing was 3.0 degrees longitude and 
2.5 degrees latitude. 
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Figure 2.2 Vertical structure of the model showing seven 
levels with dense packing near the surface. All variables 
were stored at these levels except W; the vertical velocity, 
which was found at intermediate levels (not shown, see 
igure 2.8). 
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The model is based upon the finite difference forms of 
the primitive equations which are integrated in time from 
a known set of initial conditions. The governing equations 
are the two horizontal equations of motion, the hydrostatic 
equation, the continuity equation, thermodynamic equa- 


tion, the salinity equatdon, and the equation of state 





dua ODS ane Wing + Fay (ZAB) 
dt Poa cosy JA a toS y 

da - == R m dde uSiny + Fy, 5 
3. - pş) 5 
ONY O, 105 
ST = Ayv?T + RET Y Ct CTS), (2.5) 
Ra As VS + RSS + Cs(7,S), | (2.6) 
p= po l —a(T-Te) + B(S-So)], ве.) 


where the symbols are as follows: 


t time 

A longitude 

4 latitude 

È height 

= angular speed of the earth 


iLə 





44 25 p 


So 
CTS) 


Cg(T, S) 


Ru, ks 
An, hs 


Fa 


radius of the ecarrh 

acceleration of gravity 

zonal cürrenti component 

mecidiona C cürrent component 

horizontal velocity 

vertical speed 

temperature 

constant reference temperature 

salinity 

constant reference salinity 

rate of change of | due to instantaneous 
convective adjustment 

rate of change of S due to instantaneous 
convective adjustment 

density 

density of water at reference salinity So and 
temperature To 

coefficient of thermal expansion of sea water 
cociinomaonnaoi saline conuracozon Qanun 
vertical eddy diffusion coefficients for heat 
and salt, respectively 

lateral eddy diffusion coefficients for heat and 
salt, Ze Lİ 

eastward İTTCLT007 “19570 00000 0 7 
northward frictional Loree sper inde omens 


horizontal gradient operator. 
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The model circulation is determined primarily cb 
surface wind stress, and to a lesser extent by them lune ome 
heat and salt at the surface. Figure 2.3 shows the annual 
facta tludine Tame ane Zone tes Gite иси AO CE SMS a 
(Hellerman, 1967) which was used in this model. The distri- 
purion is practically symmetric except for the ceciimrerce som 
a single minimum of the easterly stress north of the equator. 

The net heat flux at the upper boundary is calculated 


according to Haney (1971) from the equation 


Q” Qış Aal Ta- Ti), (2.8) 
where: 

Q 5 16 net sheac. flux Amen tne soceau, 

Qí is that component of heat flux that would exist 


if there were no air-sea temperature difference 
(this is primarily radiative fluxes), 

Qe js «he contribution due to the existenceuer a 
gemperature dilfterentizi"”betyeen 56 06507 
the atmosphere (this is primarily sensible and 
latent heat fluxes), and 


TA is a prescribed surface air temperature. 


Haney's results, which were based on climatological data 
after London (1957), were averaged across the equator to 
obtain the symmetric patterns of Q,, Qə. and la shown in 
figure 2.4. This form of boundary condition essentially 


couples the ocean surface temperature to a prescribed 


in 
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Figure 2.3 The annual latitudinal mean zonal stress for 
all oceans combined after Hellerman (1967). 
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Figure 2.4 The atmospheric surface air temperature 
(bottom scale in degrees Celsius) and the heat flux 
components Q4 and Q4, in lyday^i and iyday-lideg-i (top scale). 
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atmospheric temperature. To accomplish 2 nev heat a tue 
zero requires that Ta be less than lą over at least some 
portions of the domain. 

The boundary condition on salinity is quite different 
from that on temperature because the downward flux of 
salinity is determined by the rate of removal of fresh waver 
wom the surface of they ocean. Rather Thanspeecify inguthe 
salinity at the surface, the downward flux of salt was speci- 
fied from the zonal mean climatological rates of precipita- 
tion and evaporation after Budyko (1963). These curves are 
shown in figure 2.5. The data were averaged across the 
Smmator to obtain the symmetric curve shown by the dotted 
mime, In this way, the only asymmetric boundary condition 
was the surface wind stress. This curve has the property 
that the meridionally integrated net rate of evaporation 
Eumus precipitation is zero and, since no flux of sait was 
allowed at the lateral boundaries or ocean floor, was suffi- 
cient to guarantee that the volume averaged salinity of the 
ocean remained constant in time. This does not imply, how- 
ever, that no variations in surface salinity patterns would 
tend to occur. Clearly, near the equator where E-P has a 
locat maximum, a local maximum of sait would tend to be 
formed. In the region of the model where precipitation 
exceeded evaporation the most, a local minimum would tend to 
be formed. Whether the horizontal and vertical gradients of 
salt caused by the upper boundary condition would modify the 
mode circulation was an important question to be answered by 


ENIS investigation.: 
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Figure 2.5 Evaporation and precipitation rates (solid 
lines) after Budyko (1963). The net rate of evaporation 
minus precipitation averaged across the equator is shown 
by the dashed line. 
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Details of the space differencing of the momentum equa- 
tions are given by Haney (1971), and only the equations of 
temperature and salinity are treated here. These are con- 


veniently written in symbolic form as 


S = SALT ADV. + SALT DIFFUSION + CONV. ADJ., (2.9) 
+ = TEMP ADV. + TEMP DIFFUSION + CONV. ADV., (2.10) 
Uere: 
SALT ADV. = = (va: (vs) + à-tw£)], (2-08) 
A 
TEMP ADV. = = İv: (VT) A |, (5500 
4 % 

SALT DIFFUSION = AsV S + hs I$ (2.13) 

1 лт 
TEMP DIFRUSION = Aw¥*T + Ru E (5 

C 


and Ag, Re, Ay, and Pu are the diffusion coefficients as 


previously listed. Other numerical models have used values 


of Дезо2 ии 


and Kn 1] maa u (Bryan, 1969; Munk, 
1950). Because such a simple form of mixing was used, it 
mes reasonable to use similar values fer these 0611167 2 
The coefficients of vertical eddy diffusion of heat and salt 
were ku, kę = .1.0 em*seo *. The coefficients of horizontal 


eddy diffusion of heat and salt were An, As аг 101 


em^sec 1. The friction terms in the equation ol movsonm WE 


22 





similarly formulated in terms of simple horizontal and 
vertical eddy diffusion of momentum with coefficients 


l and Ru = 10.0 em^sec |. This gave a 


Am= 1.0 x 10° cm sec” 
Batlo ob ky to ky of ten which appeared to give good results; 
Krivelevich (1970) found a ratio of twenty to give the best 
results in a forty layer, two-dimensional equatorial model. 
The discussion of the "CONVECTIVE ADJUSTMENT" is given 
in a later section, and only the differencing of the advec- 
imeem and diTTusion terms are treated here. Furthermore: 
since the forms of equations (2.9) and (2.10) are identical, 
and both variables § and T are located at the same grid 
points, only the finite differencing of the salinity equation 
is given. However, because the surface boundary conditions 
me different for the vertical diffusion terms of each, a 
separate treatment of them is necessary. Although all 
finite differencing was done in curvilinear coordinates, 


only the cartesian forms are shown for the purposes of 


demonstrating the techniques involved. 


B. HORIZONTAL AND VERTICAL FINITE DIFFERENCING 
qz Thel orizontal Grid 
Ae portion. ot the hori zoncalecrid of sa ze 21 salou ts 
shown in figure 2.6. The temperature | and salinity S are 
located at the V and \ erid points, whereas the velocity 
components M and v re 1003660 at intermediate points; AX 


is the E-W grid distance, and by the N-S grid distance. 


> 








Po 


isi 


vi 


Figure 2.6 The horizontal space staggered grid with 
salinity S and temperature T located at the b and 4 
grid points and the velocity components u and v at 
intermediate points; 4% is the E-W grid distance and 
au Uhe N-5S grid distance. 
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Figure 2.7 Grid points near the western wall used to 
derive the finite difference analogues along CRE 
boundary at the general grid point (1,3); 
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e.  Differencing ofarlonizonoa1x0 0 en 


Expanding the horizontal adwecitton por Tenso] 


equation (2.11) gives 
(vs) = $ (us) ae). (2.15) 


senstdering first only those erid points intersonuto the 
boundaries and using centered finite differencing at the 


general grid point (i,j) gives 


yet (2.16) 
AX hy 
Since neither u or v is explicitly defined at the half 
integer points, a linear approximation is used; for example, 


the term (S), is approximated as 


COMER E ..—... | | Stat SLI | 211) 


Bang a similar procedure for the remaining three terms in 
equation (2.16) gives the final differenced form for 


equation (2.15), 
PSJ = [aiw js, © 21100, 6-02] [ Stat + Sis] 2 (aax) 


— Һан, + Wi] | Sina + si, | / (да) 
“ts Гә + ivy, it] | Suin + sis] / (4bu) 


7 | Vita + fəndən, E E / (^b), (2.18) 


ео 





The use of a centered difference con ie 
the boundaries is not practical because the velocity compo= 
nents u and v are not specified exterior to the walls. The 
nite difference analogues are derived MSAN E 
boundary conditions and using one-sided approximations to 
the derivative when necessary. Because the mathematical 
forms of the difference analogues are similar for all four 
nns, i is sufficient to demonstrate the differencing for 
only one of them. A segment of the western boundary region, 
moe CC [pure 2. c iS Used to derive tche diller nee 
approximation at the general point (1,j). 

Since the normal component of the current u must 
vanish along the western wall, the first term of equation 


7 15) reduces to 


o 
+ (us), = COWACOM (NE 
ex M 


Further reduction leads to 


+ (us) = (мм, t UM rt | | Si) + 82.3) : (2.20) 
AA 2 AX 


To derive the second term of equation (2.15), it is assumed 


that free slip applies for the meridional speed component v, 


enat 15, 
Vo == Vy ee ane ce, UM 
Vaia = VM y yo Ҹә, | (2-227 
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In finite difference form the terms 


wi a A 
d: Au, 
Application of the boundary conditions given by equations 


(2.21) and (2.22) gives 


Ò (v$)1) = (Wie) (SU + Saja} - (Yor, i-e Siu) (2.24) 


"y "= 


| o Difrerencine of Norizontal Diftusien 
EROS Ono һә nerizonzał diitus LonN"PORCRONNGE 


equation (2.13) gives 





dis des 
As ves As | - ES Ы (2.25) 


This is rewritten as 


Asv’S = Asis 2) & [as 
- PE "hs | Ga 


Evaluation at the general interior grid point (i,j) results 


in the equation 


e 


a ә 
s vA Sls = Al lion AA 


as \ AS | 
+. ^s [S həə, » 9 | TELIZS 


by (29298 


enr 





Application of a centered derivative to:€400 137000 3 T 


gives the final result, 
Asv* Si = As| $45 728, * Sin] 
(Ax) 


ssl 331-293 5 Side 
(Ay)? 


(229) 


The differencing of equation (2.25) along the ocean 
walls is handled similar to the evaluation of the horizontal 
advection terms along the boundaries; the differencing will 
be illustrated for an arbitrary point (1,j) along the 
western boundary. Boundary conditions along the wall are 
wose of perfect insulation, 

3$ a ts | (2.29) 
Miereing equation (2.25) in finite difference form at point 


el) gives 


| 45 
As ve Saq = ME: Eu del. Ass Ala Bla] (22320) 
BA 
de 


mielicacion of the boundary conditions and further reduction 


or the first order terms gives the final result, 


ASUS ij = 2As| Sai- Sad | p As] Samen 24 ŁO] (2,31) 
(dx) (ayy 
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4, The Vertical Grid 
Figure 2.8 shows the placement of layers in the 
vertical. The vertical index k and all os theolwspa db 
25, 8, Тә e) are defined at the integer levels shown Dy 
BRE solid A whereas w is defined at the half integer 
levels indicated by dotted lines. The (KM + 1) unequal 


thicknesses are defined as 


At py, = Er tk, KE 2 cove EM) (2.32) 
KZ lh 2 -Atı, (2.330 
Qt ewe = 2 (Ekm t H), (2, 52) 


where ZŻk<O while H>0 is the ocean depth. This definition 


satisfies the equation 


KM 
nüz “ә 2 AT k-yo. da z| mas] = H. (2.35) 
2- Ka. 


The thicknesses centered at the levels at which (u, v, 5, T, 
p are located are Celanedsac (heaven Sn 000 


MEA? integer levels, that 15, 


AZK= = EZ 4 ... Kei, KM, (2536) 


) 


inus definition has the property chat 


KM 
5 һар= Н, (22218) 
Ke 


ə 











“ed 


. 
© 


WKM+Y/2= 


777777777/7/Р/?У/7777777/7” 


Figure 2.8 Location of the variables in the vertical grid. 
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5. Finite Differeneing of Venuaegi econ 
Considering the intermediate levels A232 aes әр 


the vertical advection portion of equation (2.11) is written 


at the Kah level as 


Ś (wS)k = (WS) ksh. = иб) кәм ‘ (2.38) 
ee AX 


Eunce 15 mot available at hali levels, the average of the 


salinities at adjacent levels is used to give 


2 (wsk = Wik-Ya, [Sr-1 + Sr) —- wera, (Sk + Sioa) R (21:30) 


ASA 


Eu uncupper boundary the vertical velocity is zeros 


mathematically 
WLF O. (2010) 


Equation (2.38) evaluated at K=l gives 


o 
& (wS)¢ — [B - (8) | 65.17 
42 Als 


Further expansion gives the analogue at the level nearest 


the surface as 


&(w$), = ws. [81 +63) , (2712) 
Az 


20%, 


p 





At the ocean floor Wu = O , and at the bottom level 


equation (2.38) becomes 


è tws = er Semi + Sen) (2.43) 
de | 2. AE km 
6. Differencing of Vertical Diffusion 

Considering the interior levels first, the vertical 


diffusion portion of equation (2.13) is written at the 


general Kah level as 


den Mon] 
d = Reldulkve dale], (2.44) 
dk AZ, 


Further reduction of the first order derivative terms gives 


ES ED > SK= Suet | 
Ы c Esl Aza CN y KEZ nes KML 577 
at Je bx 


The lower boundary condition of perfect insulation at 


Z=-H and index KM*V2, is 


ə = O. 000 
Jz je T 
Using this boundary condition gives 
0 
ONE 
ка = Ы aR / kunt AU krą Mo | (2.17) 
darin 


BÈ kM 


De 





and the analogue for the lowest level becomes 


Ку 225 ) m "P+ Sian | 3 (zal 
art EM Һӧкм" b E kt - Va, 


Up Ge this paint the- -finite differencing or WNE 
salinity and temperature equations has required only a 
single treatment as both equations have similar mathematical 
forms and boundary conditions. The difference analogue 
for a term in the heat equation is obtained from its counter- 
part in the salinity equation by replacing S with the 
variable T. At the ocean surface this procedure does not 
apply since the boundary conditions are very different. 

The upper boundary condition for the heat equation 


is 


kr AT | 20 
elka “RoC Fr 


where Q=Q, + Q2(Ta-T,) Equation 2200015 ei ii”? 
strong condition that ties 14, the temperature at the 

first level below the surface, rather closely to the speci- 
fied atmospheric temperature la. Should some mechanism 
cause a large difference between the two, heat fluxes are 
ereated that tend to restore chę dlfference to a small) Vane. 


The boundary condition for salt on the other hand is 


Rg AS E 0 (2.50) 
S een. SUE ?) 


S 





This is a relatively weak condition of nearly ə 
flux which tended to make the model ocean respond very 
Slowly to the surface salt flux. Areas where the rate of 
evaporation exceeded that of precipitation tended to get 
more salty until advection and/or diffusion became suffi- 
cient to create a balanced state. Similarly, areas with 
excess precipitation tended to become less salty with the 
approach toward steady state being very slow. The reasons 
for the slow salinity relaxation time and its effects upon 
the dynamics of the model are thoroughly discussed in 
section Lil. 

The finite difference form of the vertical diffusion 


term in the salinity equation at level K*à is 


CEN 
RBS) - Reel dein dd] Ae | (2.51) 


Use of the boundary condition (equation 2.50) gives 


Rss) = [ SiE-P)~ ks 8 wi (2.52) 


JR? Һә. 


which reduces to the final form 


ks te | = Si(E-P) = Rs (S1- 82) 5 (24583 
RIL bi, KE, *AŻY2 


In a similar manner the difference form of the 


vertical diffusion term in the heat equation at the level 
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Kei, is 





СӘ 
Rr ÙT Y = Rra Фел), (2.54) 
dw ’/ı A 


Using the upper boundary condition for heat (equation 2.49) 


yields 


br z Qt Qr (Tr =T ie kr [M-Tu . (2.55) 
L 


x PoG AZ A. bts 


C. EQUATION OF STATE 

|.» Density as a Function of Temperature and Salinity 

Density depends ina nonlinear fashion upon tempera- 

mare, Salinity, and pressure; however, the assumption of 
ineompressibility reduces the dependence to two variables, 
salinity and temperature. If Po is the density correspond- 
ing to arbitrary reference values of temperature le and 
Salinity So; it follows from a Eoydler 0 expansion 
about İn and Se that the new value of density at salinity 


Sot 4S and Tot&l is given by 


"n 
pl Sotas, TotaT) = p(Sa;To) + SPIS ede AT g asa L ke әта). (2.56) 
As eT ni} 95 Ar 


Data taken from Von Arx (1964) and reproduced as 
figure 2.9 shows that over the range of salinities and 
densities commonly found in equatorial oceans, the nom eae 


terms may be dropped with little sacrifice in accuracy and 
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Figure 2.9 Plot of sigma-t vs. temperature and salinity 


2 “r Von irx (1961). 





the equation written as 


e Po | L-a(T-To) + pi S- Se) |) 


where 


(4502 

ec the density at temperature T and salinity S, 

do x Woe Cleo Sty Sie rr r) SPS Oe 6 0 0 rE Te ana 

salinity So, 

E" ebrios. | 5. 

= 190 

B= m Te CONST. (2.59) 
2: 


Solutions for İd and 8 
Introducing the quantity 


Te İOsiT,o - 4) x 105 


where the subscript zero refers to sea level pressure, gives 


a L dör 53 
u fo sr /S* MAS ) (2.60) 
B ES -— 
Po a 


to” 
Ta CONST ^ 


CZWORE) 
w 
Both A and f are determined from figure 2.9 by graphically 


taking a finite difference approximation to equations (2.60) 
and (2.61) to give 


gda = I Em x IS” u -L (Se 68), ic? 
Po wT Jseconst 


= iE 
Po (Te Ty = 2,1910 dez ) (2762) 





be rr) X M no L (E i "Sexe (sy (2.63) 
To CONS (Ta -Tr) 

THREE SS O A cael Ss ane 102,76 
em-em >, ә and 34.72%, respectively; this reference 
state is marked by a star on figure 2.9. The approximations 
Bor A and Y are taken, however, at a location on the T,S 
diagram which is representative of the temperature and 
salinity found near the surface in equatorial oceans and 
che dependence of density on these variables is nearly 


linear. 


D. INSTANTANEOUS CONVECTIVE ADJUSTMENT OF DENSITY 

Whenever a heavier water mass lies above one of lesser 
density, an unstable situation exists and natural buoyancy 
forces cause vertical mixing and a return to a neutral 
55.6. Since from the equation of state, density is a 
IENSCLIODn of temperature and salinity, any mechanism that 
eeuses temperature to increase or salinity to decrease with 
mewn is a potential source of instability. In this model 
pne only heat flux into the ocean occurs at the surface. In 
a steady state, the surface averaged heat flux must be zero. 
Thus there are some regions where the net heat flux at the 
surface is downward and some regions where it is upward, 
and it is in these regions of Upward” heat flux thay Therma, = 


induced instability is found. Salinitey-induced= iis vaste, 
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occurs in regions of the model where evaporator rnme CE 
precipitation causing a 1ayer 017060 0 006 7 7 
exist near the surface. The vertical juxtaposition due to 
advection of water masses of equal temperature but differ- 
ent salt concent is also a potential source of salinity- 
induced instability. 

Instability in the real ocean lasts only a short period 
of time. Buoyancy forces produce large vertical velocities 
which mix the unstable layers very rapidly and efficiently, 
and quickly establish a uniform vertical distribution of 
mixed quantities. The following is a reasonable method of 
modeling this phenomena. 

Figure 2.10 shows a vertical section of the ocean where 
the density increases with depth and represents the stable 
case where no convective instability will exist. Using 
previously introduced notation, several levels Kei, K, and 
K*L and the vertical thickness AÉK are shown. The inter- 
mediate levels K-3/;,, K-V}, Kev and K*9/ are shown as 
dashed lines and thicknesses corresponding to Ke V» and 
K+'/o are denoted as AEk-A4, and AŻk+'/ą respectively. 
Although in the real ocean the profile can be nonlinear, a 
linear profile ås denoted by the solid lines in figure 2.11 
Æ assumed in The model; thi” first order approximation is 
consistent with others used elsewhere and is sufficient to 
guarantee conservation ot the quadratico rity m e 


vertical advection terms. Additionally, it is assumed that 


By 
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Figure 2.11 Model Equivalent of the Staple Density PRONINE 
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the density which is storedvat each 0 ии 7 - 
vertically integrated average of the entire b e layer. 
1. фк-д ә фк for any two consecutive levels, then 
convective instability on a scale which is too small to 
be resolved by the grid will occur. Since from equation 
7 51) density is a function of both salinity and tempera- 
ture, either one or both may be the cause of the instability. 
Figure 2.12 shows the case where density instability 
is due to temperature; salinity is assumed to be constant 
with depth. Trey represents the temperature for the upper 
layer of thickness Atki and Tk the temperature for the 
lower layer of thickness AtK. To insure stability, new 
temperatures Dent and TE are determined such that a 
Slightly stable temperature lapse rate results. Mathemati- 


cally the following must hold: 


* JA 
Thea — In ifs 


* 
Sy ar Tk- e Tk + Er A £ k-və. (2.64) 
NE koVo_ 


where €T is a small positive constant whose value is the 
new lapse rate shown as the short dotted line in figure 
2.12. To insure conservation of heat, the thickness 
weighted average of the temperatures before and after mixing 


must be the same, that is, 


Tus Ara + TRAD = un BEL t Tes azn (2.65) 
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Figure 2.12 Diagram of an unstable state due to the 
temperature profile showing the necessary temperature 
adjustments to achieve stability between levels Z and 
ҸЖк-1. The solid line is before adjustment, the dashed 
EB alter. 


42 





Equations (2.64) and (2.65) comprise a set of two simultane- 


ous equations in two unknowns and when solved give 


m. s TRAEK + l Ween  GGAE k- V2.) AZ ket ) 
Ak- + A?k 


5 


* 
peso t. Tk + Er AGRA, (2.67) 


Equations (2.66) and (2.67) define the new temperature values. 
If the temperature is constant but a decrease in salinity 
vith depth causes instability, convective adjustment must 
take place such that the salinity values $& and $«., are 
replaced by new ones ar and Cur such that a stable, 
Abt ly positive vertical salinity gradient, Es, is 
established. Assuming conservation of salinity, and using 
527 .--oningsimilar to that used for temperature, gives the 


following results analogous to (2.66) and (2.67): 


SM A Św lx + ( Śk-1 + EsAŻk-V2) AS Key | (2.68) 


Һе, + AŻ 


Ska m Sr = Echte. (2.69) 


For an incompressible ocean, the value for & is zero, 
but it is reasonable to define it to be some small positive 
pmalue to avoid possible instability due to round off error 
in the computation; €r is defined in terms of To, the 


surface reference temperature, and H, the depth of the model 


ocean, as 
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era OT, (2.70) 
H 
Having fixed €x, the appropriate value for er ondes 
can be derived from the equation of state and the definitions 
of Es (2.64) and €r(2.70). If it is assumed that the 
salinity and temperature have equal weights in determining 
the density lapse rate, then it can be shown from the equa- 


muon of state that 


Es = ra ел, Cae 


A very brief explanation of how the convective adjustment 
is programmed is as follows: 
(1) The density lapse rate, (Urna — ex) / aze-v.) 55 
computed for a layer and compared with the stable value ©. 
(2) If the density lapse rate is more stable than Ер, no 
mange 15 made in temperature or salinity; if it is less 
stable than €p, the temperature and salinity are both 
"Oiusted using equations (2.66), (2.67), (2.66), and (2.69). 
(3) All layers are processed in order from the first below 
the surface to the last above the lower boundary. In the 
Convective adjustment of a pair of levels 1571007005 ИРИ. 
to render an adjacent pair of levels unstabie; therefore, 


the adjustment proceeds in repeated scans until no 


LH 





adjustment of temperature or saliniey 15 pertermedzaa sean 
scan (see figure 2.12). The convective adjustment is then 


terminated. 
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Peels TIME DEPENDENT BEBAVIOK Or Sie ИЛИ и 
1. Salinity Relaxation Time 

The upper boundary condition for the salinity equa- 
tion was the flux of salt into an initially isohaline ocean. 
However, some insight into the time dependent behavior of 
the model can be gained from the analysis of a more idealized 
кии јол, опе where the upper layer of an initially salt- 
free ocean of infinite depth is raised to a constant value 
s” Neglecting the effects of advection, horizontal diffu- 
sion, and convective adjustment the salinity equation (2.9) 


and its initial and boundary conditions become 


as — Rs XS 150 
dt Jue 

S (4,0) = CON STADT = gó Sor t70, (5202) 
S(o,2)* O. (3.3) 


To solve equation (3.1) two new variables T and f are 


defined as 


ә AE SĄ (3.5) 
NT 24 ast 
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Transforming (3.1) into T and q space gives 


z —_a (3. 


Assuming that S depends on time t only through mg implies 


that 


AS dS 46 _ 48 - È 





75...” 5”, — — Ge 
IL OG Ot = | vem 
or that 
us -— E a (3. 
e at \ er 
Substitution of this result into equation (3.6) gives 
PS I - Zr as © (3 
eg er 
Letting qa 99 yields 
T 
$$ — — Ce) =. 
E Y 
whose solution is well known as 
— 
= фә : . (5 
Further integration of (3.9) gives 
| " a 
x –Ҹ̧ 
S= ge f qe dy (3. 
O 


7 


6) 


7) 


8) 


9) 


10) 


11) 


YA 





where 
St= S(gs0)- S(3*0), X»0. (3.13) 


ranally the solution becomes 
© Se 
Sic) = sk, $ qe de (3.14) 


* 
where S is the value of § at (50. 


The value of Qo is determined from equation (343) аб follows: 
St a) x St) O for 870, (9515) 


ma@om equation (3.14) 


Sla) = CREER (3.16) 
since fora > /=7/2, (3.17) 
the final solution Гог 4 becomes 

кк == 5" (3.18) 


Jw 


and the solution of equation (3.1) including the boundary 


condition becomes 


Sig) = S*(1- ERE(C)) (3.19) 
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where 


ERR) = £ (enag, (50 
İm 6 
Higune 3- llis aplot of eguatron (3: TOE wa 

used to determine an estimate for the salinity relaxation 
time. From this graph it was clear that Кәбә gives Se $ 
the true steady state solution; however, since both the 
ECHO ME . and the vertical diffusion coefficient for salt, 
Rs, were nonzero, a true steady state would exist only as 
mume approached infinity. To avoid ambiguity, it was 
convenient to define a quasi-steady state to exist when 
S v0.9 S^ which corresponded to a g value of 6.1. Solution 
of equation (3.5) for % gives the time required for the 


level @ to be affected, 


t= 27, (3.21) 

Akor? 

Since t varied directly as the square of the depth, 
the salinity relaxation time was expected to be very long 
for the lower layers. For Kg = 0.1 emósec”* and = 90.1 
the first layer at five meters computed to be steady state 
mm 0.2 years, the third in 20 years, the riftn in 5007728157 
and the seventh in 32,000 years. As the model was inte- 
grated over the relatively short period of 30 years, theo- 


retically only the first 90 meters Were expec vcd vo bem 


equilibrium, Furthermore, sinec ume upper boundary condit ron 
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Figure 3.1 A plot of S$ vs. Ga R/(2VRE); the dotte 
line shows a quasi-steady state value of O.ĄS% 
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vas the specification of the flux of salt andi nor th: 
idealized condition treated above, the relaxation time was 
expected to be even larger for each level than indicated 
Bu (3.21). 

Figure 3.2 is the salinity anomaly vs. time plot 
at various levels for a 30-year period at the grid point 
I1=lo, J=U; (see figure 3.20). This location was chosen 
1507 studying the dependence ol the salinity relaxation time 
upon the vertical diffusion because the prescribed upward 
ix Of salt at the surface was a maximum here. düddition- 
ally, because both the horizontal and vertical currents 
were small with little shear, the effects of advection were 
expected to be negligible. The apparent breaks at approxi- 
mately 10 years in the otherwise smooth curves were caused 
by a change in Ru; the vertical diffusion coefficient or 
momentum, and did not affect the general nature of the 
adjustment. As predicted from the theoretical analysis, 
levels above 100 meters appeared to be approaching steady 
state and showed little change during the last 10 years of 
integration; however, levels below 250 meters showed a 
constant rate of increase of approximately 0.002'he per 
century at the end of the integration. Although the salinity 
at the lower levels was not in steady state, it was felt that 
this was not a serious problem because the salinity gradient 
did not have a significant effect upon the model dynamics, 
which were primarily controlled by the surface wind stress 


and the near surface ecdy momentum transie preocecscar 
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257 rc 2.24 A plot of ‘salinity vse time at the geridrpoint 
with maximum upward surface flux for levels 5, 20, 100, 
250, and 700 meters. The ordinate values have been shifted 


wo place all curves on one graph and only the siopes are 
5 enificant. 
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2. Tempenatvuremre |axat ion Time 


As discussed in Section II, the upper boundary 
condition for temperature was quite different than the one 
for salinity. The flux of salt at the ocean surface was 
specified as a function of the existing salinity and the 
difference between the fixed rates of evaporation and pre- 
Би ication (equation 2.50). If there were no horizontal 
2 sport of salt, a steady state could not be reached. 
Because the prescribed flux is nearly independent of salin- 
ity, those regions with precipitation excesses would tend 
pembecome less salty, whereas those with net deficits would 
BRASCO become more salty. Temperature, on the other hand; 
was different because the ocean surface temperature was 
tied to the apparent atmospheric equilibrium temperature 
TA TA +O4/Q2 (equation 2.49). If the temperature of 
the near surface layer Ty became much different than Ta% 
then a large heat flux resulted which drove Tą towards TA" 
resulted in a reduction of the heat flux. If there were 
no horizontal transport of heat, a steady state could even- 
tually exist in which the temperature of the ocean was every- 
where equal to age For these reasons, it was apparent 
that the theoretical analysis just performed applied more 
closely to temperature than to salinity. Accordingly, it 
Was Expected that the time requiredito reach temperature 
equilibrium would be less than that needed for salinity to 


reach steady state. To decrease the temperature relaxation 
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time, an initial temperature lapse Tat or -1°/100 meters 
was chosen, a value that was approximately equal to the 
steady state lapse rate in the equatorial domain of Haney's 
(1971) model. 

Figure 3.3 shows the temperature vs. time plots for 
several levels at the same point as in figure 3.2. As 
expected, the near surface levels reached steady state | 
comparatively quickly, whereas the deeper ones did not reach 
equilibrium over the 30-year integration period. The upper 
four levels achieved equilibrium after only 10 years and 
then showed little change for the next 20 years, whereas 
the corresponding salinity values, although approximately 
sweady state, still continued to show some rate of change 
over the corresponding period. At the end of 30 years, the 
bottom ocean levels appeared to be warming at a constant 
> of 0.25 C per century. This same trend was noted in 
Bryan's model (1968) although his time integration was 
carried out over a much longer time span. As in Bryan's 
model, the gradual heating was probably due to the selection 
of too cold of initial temperatures for the lower layers. 
Again as for salinity, the failure to reach steady state 
in the bottom layers was not felt to be especially damaging 
to the dynamical results obtained. The circulation was 
determined to be primarily wind-driven, and it is the 
surface layers which are of major interest in an ocean 


forecasting model. 
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Hipure 3.3 Temperature vs. time plots for levels 5, 20, 100, 
250, and 700 meters over a 30-year period. The ordinate 
values have been transformed various amounts to place all 
curves on one plot and only the slopes are significant. 
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3. Dynamic Current Relaxation Time 

Munk (1950), Charney (1959), Veronis (1959), and 
other researchers have proved that the wind plays an impor- 
Eat role in determining the large scale features (oi seme 
ocean circulation. It is generally known that the negative 
zonal wind stress in the near equatorial region causes a 
strong surface current towards the west, and that Ekman 
divergence is responsible for a meridional component of 
circulation at the surface towards each pole. Additionally, 
it has been established that the ocean response time to 
the surface stresses is very short, approximately 100 to 
200 days. 

No formal mathematical analysis of the current 
pesponse time is attempted here, and only the actual results 
of the time integration are presented to verify that the 
dynamic response time is small. Figures 3.4, 3.5, and 3.6 
are plots of the zonal current taken at the equator midway 
between the eastern and western boundaries. The computer 
гед data was taken from the time integration every 25 
days. Figure 3.1 is representative of the second level 
below the surface where a strong countercurrent toward the 
east of 180 cm/sec had developed during the first 10 years. 
Figure 3.5 shows the flow at 100 meters where a sizeable 
countercurrent continued to exist after adjustment of Ry. 
Figure 3.6 is the zonal current vs. time at 700 meters and 


represents typical deep ocean water. The response time of 
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Figure 3.4 Apot of Zonal current vs. time at the gaia 
point I=10, Ta Ą located on the equator midway between 
the eastern and western boundaries. The depth is 20 
meters. 
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Figure 3.5 Zonal current mM (cm/sec) vs. time for 
grid location X*€!0, 7*9 at 100 meters. 
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Figure 3.6 Zonal current u (cm/sec) vs. time at grid 
point T#!0, J=24 at 700 meters. 
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about 100 days was extremely fast compared to those of 
salinity and temperature and appeared to be independent of 
depth. 

The rapidity with which the wind-driven circulation 
adjusts was dramatically demonstrated after hw was altered 
A correct two unrealistic features of the circulation that 
ue appeared during the first 10 years, An unreasonably 
mm@one equatorial surface current had developed, and its 
associated countercurrent was located only 20 meters below 
sea level. Within 100 days after increasing Rm, the currents 
at most levels had adjusted to new values radically differ- 
ent from the previous values. At 20 meters the flow shifted 
from 60 cm/sec from the west to 47 cm/sec from the east. 
naGitionally, the core ol the countercurrent now appeared 
at 100 meters, a depth more representative of its position 


in the real ocean. 


B. EQUATORIAL UNDERCURRENT 

One of the unique features of the Pacific Ocean is the 
existence of a thin, swift countercurrent which is symmetric 
about the equator and moves eastward below the westward 
flowing South Equatorial Current. This phenomena, called 
the Cromwell Current, was first studied in detail by Knauss 
(1960) who determined the current to be about 300 km wide, 
approximately 200 m thick, and with velocities as high as 
150 cm/sec. The core, a jet-like structure with stronger 


shear above its center than below 1t; was found to extend 
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over at EEC 3500 n mi at a depth that varicd fron moo 
meters at 140 W, to about 75 meters near the Galapagos 
Islands. In a later study Knauss (1966) found a much weaker 
flow, and concluded that the countercurrent probably exhib- 
ited considerable year-to-year variability. He theorized 
that as a first approximation, ehe countercurrent scoutdarc 
considered as a separate flow imbedded in the surface South 
Equatorial Current. A stronger than average surface current 
was thought to be accompanied by a weaker undercurrent, 
meereas a weaker surface current was believed to be 
associated with a stronger than average undercurrent. 
several theoretical analyses of the Cromwell Current 
have been performed. Charney (1959) showed that for a 
EF nogeneous ocean driven by a uniform wind, the zonal current 
Profile was parabolic with the maximum velocity dependent 
upon the coefficient of vertical eddy viscosity M. For 
large values of v the easterly surface current extended to 
the hypothetical level of no motion. As V was allowed to 
1: .“5c, the countercurrent began to develop; however, at 
very small values of Y , the countercurrent became so well 
developed that the flow at the surface had reversed and 
become westerly. Veronis (1959) offered a simplified 
explanation of the undercurrent using the idealized basin 
Shown as figure 3.7. He reasoned that a surface wind blowing 
from the east caused an accumulation of water at the western 


boundary, which caused the pressure gradient to drive the 
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Figure 3.7 (After Veronis) A cellular circulation 
develops when a nonrotating basin is subjected to 

a surface wind stress; the increased pressure at-the 
western boundary causes an eastward countercurrent 
to exist below the westward surface current. 
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subsurface water towards the east. Arthur (1959), by 
neglecting all terms in the x-component equation at the 
Eguator except the pressure gradient and frictional stress 
"Nes, showed that a representative zonal) profile containing 
75 ntercurrent could be calculated from values of Sur face 
wr ent velocity, wind stress, and pressure gradient ifor a 
constant w. Additionally, he concluded that although the 
solution was stable, the nonlinear terms dia and wda vere 
baby not negligible. wWUntortunately, ine lusionyonge he 
nonlinear terms would have made the analytical analysis 
“50Tficult. 

İn an equatorial prediction model vith good vertical 
5osolution near the surface, the equatorial countercurrent 
should be readily resolved. Figure 3.8 is an E-W cross 
section of zonal current taken along the model equator at 
the end of the 30 years of time integration. A counter- 
current of about 20 cm/sec had developed at a mean depth of 
about 100 meters. This current appeared to be slightly more 
intense and at a greater depth near the western boundary 
than near the eastern boundary, and exhibited considerably 
more vertical current shear above the core axis than below 
it. Additionally, at approximately 700 meters the flow was 
toward the west, as it was at the surface. For comparison, 
the corresponding equatorial cross section for the Pacific 
after Knauss (1960) is shown as figure 3.9. Although the 


observed data exhibited considerable station to station 
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variability, all of the gross features are easily identified 
in the model. However, significant differences between the 
observational and model results were found to exist. The 
velocity of the Cromwell Current averaged about 140 cm/sec 
compared to the model mean value of 20 cm/sec. Also, the 
mean upward slope of the current core toward the east was 
somewhat greater in the Pacific cross section; this may be 
BIE Lo insufficient vertical resolution in the model. 

AS a further comparison, the vertical profile of the 
east-west averaged current in the model was plotted, (figure 
3.10), with a’ typical Pacific profile determined by Knauss 
(1966). Although the magnitudes were quite different, the 
£ “al current shear above the core of the undercurrent 
was quite similar, and both were approximately parabolic 
in agreement with Charney's theoretical results. Below the 
point of maximum velocity, the real data retained its para- 
bolic shape, whereas the mathematical form of the model 
profile is not easily determined. One possible explanation 
for the difference in the profiles is that the mean surface 
winds during the observation period were below average. 
Then according to Knauss's concept of the undercurrent as 
a flow embedded in the surface current, the weaker-than- 
normal surface winds would have caused a below average 
surface flow and an above average countercurrent. Although 
this simple argument has some merit, a more likely explana- 
tion lies in the way vertical diffusion of momentum was 


handled in the model. The assumption of a constant kinematic 
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Figure 3.10 Plot comparing the average model zonal current 
shown by the dashed line with the smoothed zonal current at 
118°W in the eastern tropical Pacific as determined by 
Knauss (1966). The model time is 30 years. 
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eddy viscosity, Ru, implied that the eddies that transfer 
momentum vertically are independent of depth. In the near 
surface layers of the ocean, Rm may increase toward the 
surface due to increased wind mixing, but the exact func- 
tional dependence of Rw on È is not clear. Also, Rm should 
depend upon the vertical stability of the water column as 
measured perhaps by a large scale Richardson number, but 
again the best formulation of the dependence is not known. 
The jet-like structure of the model undercurrent was 
well displayed in a N-S cross section of zonal current 
taken approximately mana between the west and east bound- 
aries (figure 3.11). The strong wind-generated surface 
eurrent, its associated undercurrent, and the pronounced 
vertical current shear were all readily apparent. These 
Ssmonificant features were also apparent in figure 3.12, 
“a typical cross section taken by Knauss on his 1958 expedi- 
tion, during which the current was the stronger of the two 
expeditions. Also shown on figure 3.12 is the thermal 
Buructure of the undercurrent. Although the isotherms 
showed considerable small scale spatial variations, the large 
scale pattern revealed that the isopleths were concave 
downward above the current core, whereas below it the con- 
cavity was reversed. In an attempt to identify a similar 
thermal structure in the model, the model countercurrent 
was superimposed on the N-S cross section of temperature 
(figure 3.13). A cold region above the undercurrent axis 


was evidenced by isotherms that were concave downward. 
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However, this same pattern persisted for 200 meters below 
the axis, where only a slight reversal in curvature was 
noted. The failure of the model to adequately resolve this 
thermal structure is most likely due to the simple form of 
the kinematic eddy viscosity. The thermal structure in 

the vicinity of the equator in figure 3.13 is adequately 
explained by vertical advection by a strong upward vertical 
velocity above about 300 m, and weaker downward motion below 
300 m (figure 3.23). The equatorial isotherm spreading in 
the real ocean (figure 3.12) is probably caused by enhanced 
vertical mixing in the strong shear regions above and below 
the jet core. An important goal of any future research on 
this model should be to determine a more adequate vertical 
eddy diffusion coefficient, RM ; one that is a function of 


Æ tical stability, vertical current shear, and depth. 


C. CURRENT CIRCULATIONS 

The model was driven primarily by the prescribed surface 
wind stress Tx shown in figure 250 The existence of the 
single minimum of Ty displaced slightly north of the equator 
important, for it causes the curl of the wind stress to 
be greater in the northern domain. As Munk (1959) showed, 
this asymmetry would tend to cause a stronger circulation 
in the Northern Hemisphere, and the non-zero wind stress 
curl on the equator would cause a circulation gyre to exist 
there. To more fully discuss this last point the ZENEBEM 


circulation of the model as a whole will be described. 
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Figures 3.14 and 3.15 are the zonal and meridional 
current components at five meters and illustrate the general 
flow patterns. Near the equator a strong zonal current 
low towards the west had been created. Because the coriolis 
force is exactly zero on the equator, there was little tend- 
ency for this current to flow laterally away from the 
equator and the meridional current component was essentially 
zero. Several degrees on each side of the equator a geo- 
eaeeopnic plus friction balance tended to take place and 
Bun divergence was evident. The zonal component decreased 
imn magnitude and the meridional current increased as the 
surface mass transport tended to become essentially poleward 
in both northern and southern hemispheres. In response to 
the divergence of surface water, vertical upwelling of cool 
bottom water took place in the equatorial region to maintain 
mass continuity. The weak southward directed flow on the 
equator in the eastern half of the basin, and alternating 
north-south flow in the western half, is a part of the 
horizontal gyre driven by the curl of the wind stresses 
mentioned above. 

Figure 3.16, a horizontal cross section of vertical 
velocity, w, showed the existence of a strong band of upwell- 
ing centered on the equator. This figure also showed the 
effects of the boundaries on the flow pattern. Along the 
eastern boundary, strong upwelling was necessary to replace 


the surface water carried west by the strong surface current; 
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near the western boundary the large accumulation of surface 
water due to the surface current caused large downwelling. 
Near both northern and southern boundaries, strong down- 
welling was necessary to compensate for the poleward trans- 
port of surface water in both hemispheres. The nonuniform 
255:cture of the horizontal cross section of meridional 
current (figure 3.15) is caused by the presence of Rossby 
waves in the vertical mean part of the current that have 
propagated towards the west, impinged upon the western 


boundary, and become standing waves (Munk, 1950). 


D. THERMAL STRUCTURE 

Since the model was driven primarily by the surface 
wind stress, it was anticipated that the patterns of current 
flow would have a lares effect upon the thermal ScGRCLURER 
BB rc 3.1f, a horizontal cross section of temperature at 
five meters, shoved this to be the case. A relatively cold 
tongue caused by upwelling was formed at the equator. This 
minimum was most pronounced at ne en boundary where 
cold water was advected to the surface by the strongest 
vertical velocities. As this cooler water was transported 
to the west, it gained heat and the thermal gradient weak- 
ened. Two warm tongues were formed at 6 N and 9 S due to 
surface heating and weak sinking of the cool equatorial 
water after it had diverged from the equator. Further pole- 
ward movement of this equatorial water caused cooling as 


the amount of surface heating began to decrease. 
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A realistic modeling of the vertical thermal structure 
is essential for any ocean model that is to be coupled to 
an atmospheric model for long range prediction purposes. 
Figure 3.18 is an E-W vertical cross section of temperature 
along the equator and showed a thermal gradient that was 
gravitationally stable. Except for the layers above 100 
meters, which exhibited a nearly neutral temperature gradi- 
ent, all layers showed a large stable lapse rate. The ini- 
tial temperature gradient of AK per 100 meters had not 
changed appreciably over the 30 year time integration for levels 
below 200 meters. What appeared to be the top of the perma- 
nent thermocline existed at depths that varied from almost 
zero near the eastern boundary to about 100 meters near the 
western wall. The near isothermal temperatures normally 
found in real oceans at great depths had not developed 
appreciably. However, towards the end of the integration 
the two lower levels were warming at the rate of 0.25 °C per 
5: 6, and it vas felt that over a very long time period 
Mae sired result could have been realized, 

To further examine the thermal structure, vertical pro- 
files of temperature vs. depth at the three selected stations 
marked on figure 3.18 are shown as figure 3.19. The first 
station at 10° west of the eastern boundary showed no iso- 
thermal layer; equatorial upwelling had tended to force the 
characteristics of the deep water on the surface layers. 

The next station located at the midpoint between the bound- 


aries showed the development of a near isothermal layer due 
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to the shallower vertical velocity pattern toward the west. 
Stationnumber three located 10? east of the western boundary 
showed the deepest isothermal layer. All three profiles had 
a significant gradient between layers one and two, a distance 
of only fifteen meters. In the real ocean a weak wind would 
have been sufficient to mix this layer and establish iso- 
thermal conditions. One of the weak points of the present 
model is felt to be the lack of an effective modeling of 
were wind mixing and internal vertical mixing in the 


stable case. 


ENESSALINITY STRUCTURE 

The salinity distribution was largely determined by the 
meer boundary condition for salt which depended upon the 
net rates of evaporation and precipitation shown by figure 
27 Figure 3.20 is the horizontal cross section of salinity 
anomaly, (S%,- 34.0%.) X 100, at a depth of five meters. 
Although the net rate of evaporation minus precipitation 
curve was symmetric about the equator, the surface salinity 
гл was asymmetric, which showed the strong role that 
advection and dynamics played. As shown by figure 2.5, two 
minima of evaporation rate minus precipitation rate were 
located at about 6 N and 6 S. If the ocean were completely 
motionless, two minima of salinity anomaly would tend to 
Result at these locations, as the net rate of precipivarton 
excess would tend to dilute the surface water. As seen from 
figure 3.20, both minima were shifted poleward by the 


meridional current components of both hemispheres. Between 
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these two minima a relative maximum centered at approximately 
1 S formed due to the relative maximum in the prescribed 
downward flux of salt and a contribution due Lo ver ии 
advection at the eastern boundary. The salinity minimum of 
the Northern Hemisphere was more intense than its southern 
hemispheric counterpart, a feature that must be related to 
the asymmetric surface stress Tx (see figure 2.3) as all 
other boundary conditions were symmetric about the equator. 

Figure 3.21, the equatorial E-W salinity anomaly cross 
section, further showed the effect of salinity advection. 
Since the prescribed surface flux of salt in the equatorial 
ons is upward, salinity must generally increase with 
depth. Thus, the top 300 meters are generally less saline 
mün the deeper layers. Within this upper region of lower 
salinity water there is a shallow region between 50 m and 
imeem of even less saline water due to the horizontal Convers 
gence of the meridional flux of fresh water from north and 
south of the equator (see figure 3.23). The existence of a 
ferme of minimum salinity above the axis of the countercurrent 
does not agree with the results of Krivelevich (1970) in his 
40 layer model. He found a core of minimum salinity in the 
near surface region and reported that this agreed with 


observational data. 


ERY DENSETY STRUCTURE 
Figure 3.22 is an equatorial E-W cross section of density 
anomaly (e - I.0)99« 1000. Tus similar ии. 3168 


the corresponding cross section for temperature, and its 
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dissimilarity to figure 3.21, especially Inm heU nee 
meters, suggested that temperature primarily determined the 
density structure. This same conclusion was reached by 
Krivelevich (1970) who used a nonlinear equation of state. 
Because of the strong dependence of density upon temperature 
and since the ocean was heated from above which tended to 
insure vertical stability, it was concluded that the convec- 
tive adjustment process played only a minor role in deter- 
mining the density structure. Convective adjustment was 
thought to be important only along the boundaries where 


significant vertical velocities existed. 


G. MASS TRANSPORT 

As mentioned previously, if the currents were nearly 
27 OStrophic in the vicinity of the equator, then the mass 
transport across the equator would be zero at each level. 
Figure 3.23, the transport streamfunetion, showed however, 
there was transport from the Northern Hemisphere to the 
Southern Hemisphere at the surface with an opposite trans- 
port below. The dominant flow was a two cell circulation 
in which surface water was transported poleward in each 
memisphere where it struck the boundaries, was forced to 


descend, and returned as deep water to the equator. 
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IV. CONCLUSIONS 


The model is capable of simulating many features of the 
observed equatorial ocean flow such as vertical upwelling at 
the equator, Ekman divergence and surface mass transport 
towards each pole, and the existence of a subsurface counter- 
current with a realistic meridional profile. The failure of 
the model to accurately resolve the countercurrent profile 
below 100 m was probably caused by using a vertical eddy 
diffusion coefficient for momentum, Rw; that was independent 
Bİ 5: .:cal ctability, vertical current shear, and depth. 

A better modeling of vertical mixing by using a variable km 
will probably give improved results. 

The salinity pattern was shown to be strongly influenced 
by the model currents which were primarily wind driven, 
while the salinity distribution itself appeared to have 
little, if any, effect on the general current circulations. 

mes failure of 5 modelito reaehskempepoaturcgopM cp 
steady state for the bottom layers over 30 years of time 
integration was not considered serious for an ocean fore- 
casting model, where most of the variation in parameters 
tends to occur in the near surface layers. However, the 
slow thermal response time probably helped to prevent the 
formation of the deep isothermal layer and some means of 
accurately specifying the initial temperature profile is 


clearly needed. 
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